T c suppression and resistivity in cuprates with out of plane defects 
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Recent experiments introducing controlled disorder into optimally doped cuprate superconductors 
by both electron irradiation and chemical substitution have found unusual behavior in the rate of 
suppression of the critical temperature T c vs. increase in residual resistivity. We show here that 
the unexpected discovery that the rate of T c suppression vs. resistivity is stronger for out-of-plane 
than for in-plane impurities may be explained by consistent calculation of both T c and resistivity if 
the potential scattering is assumed to be nearly forward in nature. For realistic models of impurity 
potentials, we further show that significant deviations from the universal Abrikosov-Gor'kov T c 
suppression behavior may be expected for out of plane impurities. 



PACS numbers: 74.72.-h,74.25. Jb, 74.20.Fg 

The destruction of superconductivity by disorder has 
been traditionally used to probe the nature of the su- 
perconducting state. In classic superconductors, pair- 
breaking is caused only by magnetic impurities 1 , and the 
functional form of the T c suppression, when plotted vs. 
impurity concentration or change in normal state resis- 
tivity, is known to follow the universal curve predicted 
by Abrikosov and Gor'kov (AG)-. In unconventional su- 
perconductors, ordinary nonmagnetic impurities are also 
expected to break pairs, and in the simplest approxima- 
tion where the impurities are treated as point-like (delta- 
function) potential scatterers, the T c suppression also fol- 
lows the AG form. 

As in so many other respects, the experimental situa- 
tion in the cuprates agrees qualitatively with the simplest 
notions of what should happen to d-wave superconduc- 
tors in the presence of disorder, but differs in some impor- 
tant details. For example, when Zn is substituted for Cu 
in the Cu-0 planes, T c is suppressed rapidly as expected 
for a d-wave superconductor. Here, "rapidly" means that 
the disorder-induced scattering rate required to destroy 
superconductivity is on the order of the gap scale rather 
than the Fermi energy Ep, as would be expected in an s- 
wave system. Nevertheless, the initial slope of the T c vs. 
Ap curve found in experiment is a factor of 2-3 smaller 
than the universal AG curved. This discrepancy has been 
attributed to scattering in higher angular momentum 
channels by several authors, who modelled the scatter- 
ing potential with a separable form describing scattering 
in both s- and a single higher £-wave channel*^. This is a 
simple and tractable way of including the finite range of 
the scatterers qualitatively, but is neither consistent by 
the microscopies of screened impurities in the cuprates, 
nor capable of treating the limit of extreme forward scat- 
tering, claimed to be of relevance in the cuprate case&i^. 
Furthermore, in these studies T c is calculated as a func- 
tion of the single-particle normal state scattering rate 
1/tat, as opposed to the transport rate l/rt r relevant for 



comparison to resistivity measurements. 

A further paradox was reported recently by Fujita et 
al-^, who showed that the rate of suppression dT c /dp is 
significantly higher for out-of-plane cation substituents 
in Bi 2 Sr 2 Cu0 6 (Bi-2201) and (LSCO) than for Zn in the 
same system. Since Zn is thought to be a near-unitary 
scatterer in these materials, this is somewhat mysterious 
at first sight. If one accounts for the fact that the out-of- 
plane defects are poorly screened, and may act primarily 
as forward scatterers, however, the increase of resistivity 
with defect concentration will be slowed and one may be 
able to understand this discrepancy. It is the primary 
purpose of this paper to correlate, within simple models, 
the slower rates of both T c suppression and resistivity 
increase in the case of near-forward scatterers to see if 
light can be shed on this puzzle. 

Other types of deviations from traditional AG behavior 
have been observed and require explanation. In the past 
several years, electron irradiation studies have appeared 
which are able to suppress superconductivity to zero in 
a controlled fashion in contrast to early studies which 
studied only the initial slope of the T c suppression by 
disorder 9 . Rullier-Albenque et al-£ used 2.5 MeV elec- 
tron irradiation to create defects throughout YBa2Cu307 
(YBCO-123) and YBa2Cu306.6 samples. In neither sam- 
ple was an AG-type behavior for T c observed. In the 
optimally doped sample, on the contrary, a remarkable 
linear behavior in T c vs. p was observed down to and 
including samples of vanishingly small T c . Because the 
nearly "normal" samples are expected to have a small 
superfluid stiffness, this effect was interpreted provision- 
ally in terms of a phase fluctuation model proposed by 
Emery and Kivelsonii. Deviations from AG theory at 
small superfluid stiffness are to be expected even within 
the framework of mean field theory, however, as pointed 
out by Franz et alt^, who studied the problem numeri- 
cally and included the self-consistent suppression of the 
order parameter around each impurity site. These au- 
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thors reported positive curvature tails in T c vs. impurity 
concentration rii for large rij, in contrast to the negative 
curvature in the AG plot. 



I. TOY MODEL 

To understand physically how some of these effects 
might arise, we consider a simple model of forward scat- 
tering by impurities in which all effects can be calcu- 
lated analytically. We will assume that in-plane impu- 
rities to be described by ^-function like isotropic poten- 
tials, and out-of-plane scatterers to be described by an 
extended potential in the plane; this simply assumes that 
the screened Coulomb potential created by out-of-plane 
impurities produces a "footprint" sensed by quasiparti- 
cles moving in the CuC-2 planes. 

As proposed by Kee^, we interpolate between these 
cases by considering a weak scattering potential V^g, 

which exists only at the Fermi surface \k\ = kp, and 

cuts off unless k and k' are sufficiently close 



v Q if \4>-<j>'\ < < 

otherwise 
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The diagonal and off-diagonal self-energies in this model 
are then 



T, ((j),uj n ) = -T 



UJ n ) 
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where T — irniNgVQ, with No the density of states at 
the Fermi level. In Eqs. ([2]) and ([3]), the renormalized 
Matsubara frequencies and gap magnitudes are iuj n — 
iuj n — E and A = A +Xi, and we assume a d-wave form 
of the unrenormalized order parameter, f{<j>) — cos 20. 
Note we work in units where Boltzmann's constant ks — 
1. The critical temperature is determined as usual from 
the linearized gap equation 



A = V d N 2nT c V -L / d0'/(0') 2 ^ 

^ n 27T J UJ n 



(4) 



UJ n >0 



where the d-wave pairing interaction is given by 
(VdNoy 1 = l/21n(2e 7 w c /7rT c0 ), and T c0 is T c in the ab- 
sence of impurities. If we now take 4> c <C I , and T c < T c q 
we find to leading order Ao /Q n = Ao/w n , i.e. there is 
no pairbreaking by pure forward scattering, analogous to 
Anderson's theorem in an s-wave superconductor with 
isotropic nonmagnetic impurities. 

Since we are only interested in calculating the criti- 
cal temperature we can neglect quadratic variations of 
the order parameter and therefore perform the angular 
integrations in the definition of the self-energies. After 
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FIG. 1: The factor Q(<j> c ) as a function of the maximum scat- 
tering angle (j> c . It reflects the decrease of the initial slope of 
T c suppression as a function of resistivity if one approaches 
the forward scattering limit (<f> c — + 0). 



solving the resulting set of equations for the ratio A /0j n 
that enters Eq. (UJ), we can write 



1 = V d N 7rT ( 



2T(t) c - rsin(20 c ) 



(5) 



It is worthwhile noting at this point that an expansion 
in powers of 4> c applied to ([5]) as performed in Ref. 6 
is appropriate except when T c — » 0. Following this ap- 
proach would lead to the incorrect conclusion that the 
T c suppression is linear all the way to T c = in the for- 
ward scattering <p c — > limit. Instead, it is clear that 
the expression (|5|) may be summed exactly, leading to a 
modified AG result of the form 



In 



T, 



cO 



f 



irT, 



(6) 



where the pair breaking parameter in the forward scat- 
tering limit (4> c <C 1) is of the order of </>;?: 



r(> c - sin0 c cos</> c ) 2r^ 



-kTc 



3ttT c 



(7) 



Thus we see that the effective pairbreaking rate will de- 
crease dramatically as the scattering becomes more for- 
ward in nature. This is indeed consistent with the conclu- 
sions of earlier works^, which concluded that anisotropic 
scattering slows T c suppression for fixed impurity con- 
centration. 

It remains, however, to calculate the dependence of the 
normal state resistivity on the impurity scattering within 
the same framework if one wishes to compare directly 
to experiments where pairbreaking is measured by the 
increase in residual resistivity. For an isotropic 2D Fermi 
surface this may be written as 



2m 2 
e 2 p 2 F 



dcf> dcf>'{l - cos(0' - </>))m\V(<p, <j>')\ 2 (8) 
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FIG. 2: T c suppression as a function of resistivity calculated 
within the toy model. We find a decrease of the initial slope of 
the T c suppression as a function of resistivity if we approach 
the forward scattering limit. 



leading us to the expression 
h 2T 



e 2 E F 
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(9) 



so that T c may be expressed directly in terms of the 
impurity-induced resistivity as 
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where the factor Q{4> c ) is given as 

Q(<Pc) = (0c - sin (f> c cos <p c )/ ((f> c - sin <j> c ) (11) 

It ranges from Q(Q) = 4 in the forward scattering limit to 
Q(n) = 1 in the case of isotropic scattering and its depen- 
dence on <f> c is shown in Fig. [T] The factor Q(<fi c ) is also 
directly related to the initial slope of the T c suppression 
that can be derived from Eq. (|10p as 



(12) 



T c0 S^^ c,R \T c0 h 



We see that although the rate of T c suppression (or 
resistivity increase) is certainly much slower for forward 
scatterers than in the isotropic scattering case, the depen- 
dence on the resistivity is stronger for the former class of 
defects - within our model about a factor of four. This 
may indeed explain the results of Fujita et alj^, since 
in these experiments out of plane (more forward scat- 
tering) and in-plane (isotropic) impurities are studied in 
separate samples. In Fig|2] we show the suppression of 
T c as a function of resistivity for different values of <p c as 
given by Eq. (JTUJ). While the initial slope is indeed seen 
to increase as the scattering is made more anisotropic, 
approaching the value of 4 as <f) c — > 0, we see also that 
within the toy model there is no deviation from the form 
of the AG curve, as exhibited explicitly in Eq. ©. 



FIG. 3: T c suppression as a function of resistivity calculated 
within a tight binding model. 1/k gives the range of the 
scattering potential, where s< 1 corresponds to the forward 
scattering and k > 1 to the isotropic scattering case. The 
non-monotonic behavior of the curves for different k results 
from the competing effects of a reduction of the downward 
curvature with a simultaneous decrease of the initial slope for 
decreasing k. The critical temperature of the pure sample is 
chosen to be T c o = O.Oli. 



II. REALISTIC MODEL 

To confirm the above intuition, and get a sense of the 
range of behavior possible in real systems, we now con- 
sider a more realistic model consisting of randomly dis- 
tributed out-of-plane impurities with Yukawa potentials 
Vi = Vq exp(— KVi)/ri, where fj is the distance from a 
dopant atom to the lattice site i in the plane. 1/k gives 
the screening length of the impurity potential and the 
forward scattering case (k — > 0) as well as the isotropic 
scattering case (n — > oo) are included with this specific 
choice of the impurity potential. The Fourier components 
of the screened Coulomb potential can be written as 



V{k, %') 



\Vo\ 



k-k 1 



(13) 



We use a square lattice to mimic the copper-oxide 
plane and include up to next-nearest hopping terms in 
numerical computations. The electronic band is taken to 
be efc = —2t (cos k x + cos k y ) — At' cos k x cos k y — fi, with 
t' = — 0.3i and /i = —1.1 1 in order to mimic the known 
shape of hole-doped cuprate Fermi surfaces near optimal 
doping. The T c suppression is calculated by solving the 
gap equation self-consistently: 



A,= JLVV3L Ak ' 



A/3 ^ k > k 'Cjl 



(14) 



with an effective electron-electron pairing interaction 
Vgp — V d f(k)f(k'). The choice of the symmetry func- 
tion f(k) = 1/2 (cos fez — cos k y ) leads to a d x 2_ y 2 sym- 
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FIG. 4: T c suppression as a function of resistivity in the for- 
ward scattering limit for two different values of T c o. The 
upward curvature that is visible especially for very low values 
of k diminishes for lower values of T c o and the T c suppression 
as a function of resistivity becomes more linear. 



metry of the order parameter. The diagonal and off- 
diagonal parts of the self-energy are given by 



Z (k,co) = n t J2\v(k,k r : 
k' 
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where we can neglect quadratic variations of since we 
are only interested in the region near T c . To compare the 
suppression of T c to the increase of the normal conduct- 
ing resistivity we have to calculate the resistance within 
the same model. For the transport rate r _1 (fc) we use 
an approximation introduced by J. M. Ziman^ that ex- 
pands the quasiparticle scattering rate by an additional 
scattering-in term, 
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and that has proven to be very accurate even for highly 
anisotropic transport rates^. In the expression for the 
transport rate we have replaced the full T matrix by the 
single impurity scattering potential Vkk' ■ The conductiv- 
ity can then be calculated from 



d 2 k 
(2nf 



Tkvlj (e fe ) 



(18) 



assuming a cubic symmetry of the transport tensor. 

In Fig. [31 the suppression of T c as a function of resis- 
tivity is shown for the tight binding model discussed in 
this section and T c0 = O.Oli. As we have already seen 
for the toy model, the initial rate of T c suppression vs. 
resistivity increases with increasing range of the impurity 
potential 1/k, and seems to be a robust feature of the T c 



suppression due to forward scattering processes. However 
the particular shape of the T c suppression curve depends 
on the details of the considered model, e.g. the band 
structure, the doping strength or the functional form of 
the scattering potential. In FigHl it is shown that by 
approaching the weak coupling limit T c q <§; Ep due to a 
lowering of T c o, we find a more linear dependence of the 
critical temperature on the resistivity, reminiscent of the 
linearity found in the electron irradiation experiments. 
At present, we do not have an analytical understanding 
of the origins of this quasilinearity. 



III. CONCLUSIONS 

We have investigated the disorder-induced reduction 
of the critical temperature in high T c compounds due to 
out-of-plane disorder, assuming that the out-of-plane de- 
fects act primarily as elastic forward scatterers. These 
calculations may be important not only to understand 
experiments where this type of disorder is varied system- 
atically, but also to understand T c 's in systems which 
are intrinsically disordered by the doping process. Cal- 
culating both the normal state resistivity as well as the 
T c suppression in both a simple toy model and for more 
realistic scattering potentials and bands, we found that 
although the effect of forward scattering on both quan- 
tities is smaller than for isotropic scattering, the T c sup- 
pression as a function of resistivity is stronger for forward 
than for isotropic scatterers. In the case of the toy model, 
we showed that the suppression is AG-like, with modi- 
fied pairbreaking parameter. The basic physics of the 
more rapid initial suppression of T c vs p was confirmed 
numerically using more realistic Yukawa-type impurity 
potentials within a tight binding model. 

Our result stands in apparent contradiction to ear- 
lier work comparing T c vs p for anisotropic potentials 
scatterers^, where only s and d-wave components of the 
scattering potential were retained. These authors con- 
cluded that the greater the ratio of potentials Vd/V s , the 
weaker would be the suppression of T c vs p. While this 
model is mathematically consistent, it is physically un- 
realistic and does not yield generic results. As shown 
in the Appendix, more realistic potentials may be ex- 
panded in Fermi surface harmonics, but to obtain the 
proper (weaker) enhancement of the resistivity with dis- 
order it is important to retain the p-wave component as 
well. If this is done for a generic potential, the result of 
this paper obtains. 

This point leads us to an important conclusion regard- 
ing T c suppression experiments making use of Zn to re- 
place Cu. It has been recognized for some time that 
the T c suppression rate vs. residual resistivity for Zn is 
smaller than predicted by AG theory for pointlike scatter- 
ers in a c?-wave superconductor— by about a factor of three 
in optimally doped YBCO, and even smaller in other ma- 
terials like LSCO. We have now argued that the "conven- 
tional" explanation, based on highly anisotropic scatter- 
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FIG. 5: Comparison of experimental data on T c suppres- 
sion in optimally doped cuprates. Diamonds: Nd-LSCO 13 ; 
squares: Zn-YBCO-^; triangles: electron-irradiated YBCO 10 ; 
stars: Zn-LSCO Solid lines: toy model result for isotropic 
scatterers; dashed lines: for purely forward scatterers. Top 
(bottom) panel, E F /T c0 = 10(4). 



ing by Zn impurities, is unlikely to be correct. It seems, 
therefore, that even in the optimally doped cuprates 
the deviation must be ascribed to effects of strong cor- 
relations or strong coupling superconductivity. Within 
Eliashberg theory,— the AG pair-breaking parameter is 
renormalized T — > T/(l + A), where A is the dimensionless 
coupling, while the disorder-induced resistivity change is 
unrenormalized to leading order. The physical origin of 
this effect is not clear in the cuprate context, however, 
Kulic and Oudovenkc 5 argue that renormalizations of the 
scattering vertex due to strong interactions within the t 
-J model introduce significant suppression of the trans- 
port rate induced by impurities. While this cannot ex- 
plain the effects of forward potential scattering for out of 
plane scatterers as discussed here, since the renormaliza- 
tion in their model is independent of the anisotropy of the 
scattering, it may be part of the solution to the Zn prob- 
lem. Another perspective on the same physics to explain 
this effect may involve the low-energy spin fluctuations 
known to be induced by disorder in strongly correlated 
systems^— No theoretical work which considers both T c 
and residual p is available at this writing, however; it is 
difficult to predict a priori which quantity will be more 



strongly influenced by correlations. 

We summarize the situation comparing the current 
theory to experiments in Fig. [5] Since we are primar- 
ily interested in order-of-magnitude physics, we discuss 
only the initial slopes of T c suppression measured in var- 
ious experiments, and compare to the toy model result 
Eq. (|12[) . To obtain a fit one must make an assumption 
about the parameter Ep/T c o which enters this expres- 
sion. A reasonable choice for the high-T c materials is 
Ef/T c o — 10, and this choice puts most of the data on 
Zn-YBCO about a factor of 3 higher in slope, as found 
by previous authors. The range of T c suppression initial 
slopes within the current approach is then shown in the 
Figure in gray, ranging from the isotropic result to the ex- 
treme forward scattering result. For comparison, a set of 
curves is also given for the unphysical case Ep/T C Q = 4. 
While these parameters should not be taken too seriously 
in a quantitative sense, it seems clear that a) a hitherto 
unaccounted for physical mechanism is required to ex- 
plain the data on Zn-substituted samples; and b) the 
effect of out-of-plane scatterers within the present model 
can only account for a factor of 4 or so increase in the 
magnitude of the slope. Thus our work has explained 
qualitatively the paradoxical result that the out of plane 
scatterers generically reduce T c more quickly than in- 
plane relative to residual resistivity, there is still a quan- 
titative question remaining regarding the magnitude of 
the suppression in both cases. 

The final experimental result we have attempted to 
discuss is the fascinating linear T c vs p suppression mea- 
sured on electron-irradiated single crystals of optimally 
doped YBCO by Rullier-Albenque et ali£. The devia- 
tion from the AG pairbreaking result was attributed by 
these authors to phase fluctuations according to a model 
put forward in Ref. [Ill This is plausible for the under- 
doped sample where phase fluctuations are expected to 
be strong, and for the highly irradiated optimally doped 
sample, where the superfluid density is also small. We 
note, however, that within this scenario it is puzzling and 
must be regarded as accidental that the initial T c slope in 
the optimally doped sample, due to pairbreaking, is the 
same as the final slope before superconductivity disap- 
pears. It is therefore equally likely, in our opinion, that 
the quasilinearity found by these authors is due to the 
effects of out-of-plane defects created by the electron ir- 
radiation, as discussed here, and/or the effects of order 
parameter suppression as discussed in Ref. [iH 
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APPENDIX: FOURIER EXPANSION OF THE 
IMPURITY POTENTIAL 

In previous works the effect of forward scattering on 
the critical temperature has been studied by including 
higher angular momenta of the impurity potential when 
calculating the self-energies. In this section we show the 
link between the results of our toy model and an ap- 
proach, where the impurity potential is decomposed in 
its relevant Fourier components. Particularly we want to 
point out the importance of the p-wave component of the 
impurity potential for calculating the normal conducting 
resistivity that has been neglected in previous works and 
that is the key in understanding the stronger suppression 
of the critical temperature as a function of resistivity in 
the case where forward scattering is the dominant scat- 
tering process. 

Starting with the square of the impurity potential that 
is given in Eq. ([1]) we can expand it in a Fourier series of 
its two arguments <j> and 4>' 

oo 

\V(cb,<b')\ 2 = v 2 Q V + v 2 Q J2Vk (A.l) 

k—l 

x [cos(k(p) cos(k(j>') + sin(fc0) sin(fc0')] 



with Vq — (f> c /Tr and Vk = 2sin(/c0 c )/fc7r. Using this 
expansion to calculate the self-energies near T c we notice 
that the integration projects out only the s- and the d- 
wave part of the impurity potential, leading to 

S = 2ttTV (A.2) 

and 

v (A \ pt/ cos(2<ft) 

= 7rrU 2 ; (A. 3) 

Solving for A /uj n and performing the Matsubara fre- 
quency summation in the gap equation we find 

This result is in agreement with the results of Kulic and 
Oudovenko and their parameters T s and can now be 
directly compared to the toy model parameter cj) c : 

T s = 4r0 c , T d = 2r sin(2<£ c ) (A.5) 

To calculate the normal conducting resistivity one has 
to weight the quasiparticle scattering probability by its 
impact on the resistivity leading to Eq. ([5]) . It is obvious 
that in the case of pure forward scattering the increase 
of resistivity with impurity concentration is much slower 
than for isotropic scattering, a fact that has been taken 
into account by the factor 1 — cos(<j> — 4>') that projects 
not only the s-wave part but also the p-w&ve part of the 
impurity potential out of the Fourier expansion. The 
resistivity can then be written as 

Solving this equation for T and inserting it in the expres- 
sion for the T c suppression leads to Eq. (fT2")l where Q(<fr c ) 
can be expressed by the s-, p- and d-wave part of the im- 
purity potential as Q(cj) c ) = (2U - V 2 )/(2V Q - Vi). The 
dependence of Q(4> c ) on the maximum scattering angle 
4> c is shown in Fig. [T] It shows that the effect of im- 
purity scattering on the resistivity is drastically reduced 
in the forward scattering limit due to the p-wave char- 
acter of the scattering-in term while the effect on the 
T c reduction is not as strong, since the only anisotropic 
component which plays a role in the latter case is d-wave. 



